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One of the electromechanical transducer candidates for the electronic 
telephone is a bilamellar piezoelectric ceramic* In order to meet the 
design template for transducer response in the acoustic band, 0.3 kHz- 
8.0 kHz, a controlled resonant condition must be introduced at the upper 
end of the spectrum. 

An analytical program, consisting of three complementary parts, was 
carried out in order to understand the phenomenology of the transducer/ 
support system response to acoustic loading. The three parts are: (i) a 
simple direct variational model, used to generate parametric design in- 
formation; (ii) an exact solution with a lumped mechanical model of 
the support structure, used to evaluate the effect of using different rubber 
materials in relation to the design goal; and (Hi) a finite element modal 
survey of the system, used to determine the necessary design modifications 
and to expose deficiencies in the previous models. 

I. INTRODUCTION 

Several alternative designs are under consideration as transducer 
elements for the electronic telephone. 1 One of the leading candidates 
is a bilamellar piezoelectric ceramic plate consisting of two thin 
circular ceramic wafers that are electroded on both surfaces and ce- 
mented together. 2 The disks are joined with opposing polarity so that 
the flexural response of the assembly to applied acoustic loading re- 
sults in additive voltage output. 

The design objectives for the transducer response as a function of 
driving frequency are shown in Fig. 1 . 3 The cross-hatched areas indicate 
the template within which the response should fall. Important char- 
acteristics of this template are: (i) the electrical output rolls off below 
300 Hz in order to exclude low-frequency room noise; (ii) output is 

* It has since been decided to eliminate rubber from the design and to replace 
the bilamellar transducer by a metal/ceramic combination. 
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Fig. 1 — Microphone response design objective. 



relatively flat between 300 and 1000 Hz; (in) output increases by 
approximately 7 dB between 1000 and 3000 Hz in order to compensate 
for transmission loop losses and to improve speech recognition; and 
(iv) output rolls off rapidly at frequencies above 3000 Hz in order 
to eliminate crosstalk and other high-frequency noise. The object of this 
investigation is to determine the extent to which transducer support 
damping can be used to achieve these characteristics while maximizing 
the microphone sensitivity. 

The basic approach is to design a bilamellar structure with a funda- 
mental flexural resonance near 3000 Hz — this will guarantee flat re- 
sponse up to frequencies just below the resonance — while providing a 
support configuration which will permit shaping the response around 
the resonant peak. In addition, this shaping should include the sup- 
pression of resonant response at higher frequencies. The shaping of 
the response curve at frequencies below 300 Hz is considered to be 
a manageable problem. 

One support configuration that has been tried is shown in Fig. 2. 
The ceramic disks are mounted between soft rubber "0"-rings that 
are held in place by a relatively rigid housing. Such a design has one 
serious drawback— the lack of stability of the transducer response with 
respect to slight changes in rubber precompression. This sensitivity is 
due to the large contact area increase, and corresponding increase in 
support stiffness, as a function of slight changes in precompression. 
Large excursions in support stiffness will affect the location of the 
fundamental resonance and, thereby, distort the response. 

A more dimensionally stable configuration is shown in Fig. 3. Conical 
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Fig. 2 — "0"-ring support configuration. 



rubber washers are placed above and below the disk and have flat 
surfaces in contact with the ceramic, insuring relatively constant 
stiffness with respect to precompression. The stiffness of the washer 
is primarily dependent on the thickness and height of the cross section, 
with, a secondary dependence on the cone angle, shown as 25 degrees 
in Fig. 3. This design can be easily modified in order to achieve the 
response objectives by making suitable adjustments in these parameters. 
While the concern here is with controlled response through support 
damping, other concepts, such as the addition of acoustic elements 
to the design (acoustic mass, compliance, and resistance) or con- 
strained layer damping, could be considered. 3 Economic constraints 
are paramount in deciding the most feasible concept, however, so that 
manufacturability, material availability, and unit cost are vital 
ingredients. 
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Fig. 3— Conical washer configuration. 
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In the following sections, an analytical effort that consists of three 
mutually complementary elements is described. First, a direct varia- 
tional procedure is used to calculate, approximately, the first two 
resonant frequencies of the ceramic and its support system. The attrac- 
tion of these approximate methods is their simplicity of expression 
and the consequent ease in assessing system trade-offs. The second 
phase deals with the washer as a lumped mechanical network that 
acts at an "effective support radius" and is accounted for by a boundary 
condition on the shear resultant at this location. The exact solution 
for the forced response of the transducer, as a function of driving 
frequency, is then found. The third phase consists of a resonant fre- 
quency and modal shape survey of the complete structure by using an 
elastic, axisymmetric, dynamic finite element code. These finite element 
results provide a check on previous calculations and give more de- 
tailed information on the motion executed by the support system. 

From this description the role of analysis is seen to have several facets: 
(t) as a design guide (to identify parameters and to check initial ex- 
periments); (n) as a key to understanding the phenomenology; (Hi) to 
assess hard designs through detailed analysis; and (iv) to provide 
guidelines for future designs. 

II. RUBBER CHARACTERISTICS 

Before proceeding with the analytical details, a few comments on 
the thermomechanical properties of viscoelastic support materials are 
in order. From the transducer response template in Fig. 1, the primary 
information needed is the complex viscoelastic moduli over the fre- 
quency range 100 Hz-10,000 Hz. In addition, since the transducer 
must have stable response characteristics with respect to temperature 
changes down to about -40°C and up to about +50°C, the effect 
of temperature on these moduli must be known. 

For these analyses, the candidate polymers were assumed to be 
isotropic, thus reducing the number of moduli about which knowledge 
is required down to two (e.g., the shear and bulk moduli). Also, it was 
assumed that the materials were nearly incompressible over the fre- 
quency and temperature ranges of interest (the bulk modulus much 
larger than the shear modulus), reducing the number down to one. 4 
For example, if the complex extensional modulus is known, the complex 
shear modulus is found by dividing by three. It suffices, therefore, 
to know the storage modulus, in either extension or shear, and the 
loss tangent over the acoustic frequency range at temperatures in the 
environmental range. 
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An additional simplification is possible by assuming the polymers 
to be thermorheologically simple, so that the theory of reduced vari- 
ables 5 applies (e.g., frequency and temperature are interrelated). Thus, 
if the complex extensional modulus as a function of circular frequency 
oj at a temperature T„ is given by (i = \/— 1) 



E(u, T.) = E'(u, T.) + iE"(u, T ), 



(1) 



where E' and E" are the storage and loss moduli, respectively, then 
the extensional modulus at a temperature T is given by E(coa T , T ). 
The multiplier a T is referred to as the time-temperature shift function. 
For a typical polymer, such as polybutadiene, a plot of extensional 
modulus versus frequency (see Fig. 4) at room temperature can be used 
to generate data at other temperatures provided the shift function has 
been experimentally determined and provided that the room tempera- 
ture data extends over a sufficient frequency range. 

Characterizing each polymer to this extent is prohibitive, however, 
and the usual approach is to generate data at a fixed frequency while 
varying the temperature. Such data is shown in Fig. 5 for two polymers 
of interest: (i) a blend of 50 percent cis-4 polybutadiene and 50 percent 
styrene-butadiene rubber (called PBD/SBR) and («) a blend of 75 
percent cis-4 polybutadiene and 25 percent chlorobutyl rubber (called 
PBD/CBT). A fixed frequency of 110 Hz (data obtained with a Vibron 
Viscoelastometer marketed by Imass, Inc., Accord, Mass.) was used 
and the temperature was varied sufficiently to capture the transition 
regions of interest. Note that the loss tangent, defined by 
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Fig. 4 — Complex extensional modulus, polybutadiene, T = 20°C. 
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Fig. 5 — E' and tan 5 vs temperature. 
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tan 8 = E"/E', 



(2) 



has been shown in lieu of the imaginary component of the modulus. 

In order to construct modulus versus frequency from these data, 
the shift function for each polymer must be known. Ordinarily, a T 
would be determined experimentally from several fixed frequency runs 
and their graphical superposition. It is often convenient, however, 
to assume a form for the shift function that is found to fit a wide variety 
of polymers and is called the WLF equation: 7 



log.o a T = - Cl (T - T R )/(c 2 +T - T B ), 



(••5) 



where d and c 2 are constants and T R is a reference temperature. Common 
practice is to use c, = 8.86 and c 2 = 101.5 as the constants and a 
reference temperature in the middle of the transition region. For the 
analytical work described here, T R for the PBD/SBR was selected 
as -60°C and for the PBD/CBT was chosen as -40°C. Then, with 
the help of (3), the data of Fig. 5 was converted to the form of Fig. 4 
for specific temperatures. 

These modulus values can be approximately converted to effective 
stiffness by using simple strength of materials considerations. If the 
cone angle is neglected and the washer is assumed to be in a state 
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of plane stress, then the factor which converts modulus to stiffness 
(per radian) can be written 

2ft 

(4) 



1(1 - v'r) ' 

where r is an average radius for the washer, t is its thickness, I is the 
height, v R is the Poisson's ratio for the rubber, and the multiplier 
indicates that both washers are being taken into account. The cantilever 
frustum, which does not provide support stiffness in its unconstrained 
configuration, is neglected. If v R is assumed to be 0.5 and the dimen- 
sions of Fig. 3 are used, k = 0.8, indicating that the effective rubber 
stiffness per radian is eight-tenths of the extensional modulus. This 
conversion factor will be used in the next section in order to help 
generate parametric design information. 

III. SIMPLE VABIATIONAL SOLUTION 

As a first step in the rational design process, a procedure for estimating 
the two lowest resonant frequencies of the transducer, as a function 
of geometric and material parameters, is developed. A Rayleigh-Ritz 
procedure is used for deriving these design equations. First, a functional 
is written which represents the strain energy and kinetic energy of 
the plate and its deformable supports, less the work done by the acoustic 
loading. 8 Classical infinitesimal plate theory is used (rotatory inertia 
and shear deformation are neglected) and piezoelectric stiffening effects 
are ignored. Then 

- W 2 f P (r)h(r){w(r, t)\ 2 rdr - f p(r, t)w(r, t)rdr 

+ hrJcAw(r, ,t)} 2 - WrM.Mr. ,t)\ 2 , (5) 

where r is the radial coordinate of the circular plate, t is the time, 
and «j(r, t) is transverse deflection. The flexural stiffness, density, 
thickness, Poisson's ratio, and radius of the plate are D(r), p(r), h(r), 
v(r), and r e , respectively. The effective mass, effective stiffness, and 
effective support radius for the rubber are denoted by M , , k, , and 
r, , respectively. The circular frequency is a> and the applied acoustic 
pressure is p(r, t). 

Next, an approximate deflected shape is assumed, in terms of one 
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or more undetermined parameters, and substituted into (5). This 
shape function should satisfy the geometric boundary conditions for 
the plate (i.e., those on deflection and slope) identically, but may also 
satisfy natural boundary conditions (i.e., on shear and bending moment). 
After substitution, the spatial integration is carried out; then, the 
stationary value of the functional is found through the first variation 
and subsequent solution of simultaneous equations for the undeter- 
mined parameters. Eigenvalues are found from the homogeneous sys- 
tem. The procedure has been used for clamped and simply supported 
plates 9 and is usually found to be within a few percent of exact solutions. 
The trial function for the microphone is taken to be 



u»(r)=a + « 1 [l-(y^' 



£' 



1 +v 

5 + v 



(6) 



where a and «i are the undetermined parameters (the harmonic time 
dependence has been suppressed). This trial function has the properties 
that 



and 



w(0) = a + <*i i 
w(r c ) = a , 

M T (r e ) = 0. 



(7a) 
(7b) 

(7c) 

This implies that the generalized coordinate a represents the motion 
at the outside edge of the plate and that the generalized coordinate 
a, represents motion of the center of the plate relative to edge motion. 
The boundary condition on shear at the outside edge is not, and need 
not be, satisfied by the trial function; the boundary condition or radial 
bending moment at the outside edge, which also need not be satisfied 
by the trial function, is explicitly satisfied, as indicated by (7c). Note 
that (6) has the value w(r t ) = w, at the effective support radius. 

Carrying out the steps previously indicated yields the matrix equation 
governing the system: 

(8) 



where the stiffness matrix, [K]; the mass matrix, [M]; and the load 
vector, {F}, are given by 

32( 1 + v)(7 + v) 
3(5 + vf . 



B-? 



\w, Xw, -f- 



(9a) 
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[M] = phr; 



\ + m 
I / 7 + p 



_6 \5 + »> 



+ mw. 



+ mw) 
(3k 2 + 36»/ + 113) 



30(5 + v y 



+ mid, 



(9b) 



and 



\F\ = \pA 



[b-m^-imm 



(9c) 



respectively. The dimensionless parameters X and m are defined by 



X = 



r 2 c r,k, 
D 



m = 



r,M , 
phr'c 



(10) 



and the variables p„ and r, are the uniform acoustic pressure and the 
loading radius (0 < r, < r,). 

As an illustration of the Rayleigh-Ritz procedure, consider a trans- 
ducer composed of two PZT-5A disks, each 0.006 inch thick and 0.590 
inch in diameter. The bonding layer is assumed to have negligible 
thickness. With an in-plane extensional modulus of 6.1 X 10 11 dynes/cm 2 
and a Poisson's ratio of 0.35, the flexural stiffness for the plate is D = 
1.644 X 10° dyne-cm. The effective radius of the rubber support (the 
centerline of contact with the conical washer) is taken to be 0.708 cm, 
r c = 0.75 cm, the density for PZT-5A is 7.S g/cm 3 , and the total thick- 
ness of the plate is 0.0305 cm. 

Using these data, approximate values for the first two resonant 
frequencies can be found as a function of the stiffness ratio, X, and 
the mass ratio, m. Figure 6 shows these two resonances plotted pa- 
rametrically with respect to X and m. From this plot, the primary 
effect of the rubber effective mass is to lower both resonances (the 
second much more markedly than the first). 

The Rayleigh-Ritz results are summarized in Table I. Effective 
translational inertia is found from Ref. 10, where the effective mass 
of a rubber block bonded between two plates was shown to be slightly 
larger than one-third of the total mass. Since the total rubber volume 
per radian is 0.0312 cm 3 and r, = 0.708 cm, then 



M, = 0.0146 p„ g/cm/radian, 



(11) 



where p R is the density of the rubber (p R = 1.2 for PBD/SBR and 
1.0 for PBD/CBT). The extensional modulus values are obtained by 
fitting the WLK shifted data of Fig. 5 by collocation; 11 the effective 
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Fig. 6 — First and second resonances vs rubber mass and stiffness. 

stiffness is then computed using the conversion factor, k = 0.8, found 
previously. The frequencies shown in Table I will be seen later to 
be in excellent agreement with measured results. 

The direct variational calculations can be extended to include forced 
response and complex rubber properties. Rather than rely on (8) 
entirely, however, a more exact representation is formulated in the 
next section. 



Table I — Rayleigh-Ritz Results 
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IV. EXACT SOLUTION-LUMPED SUPPORT PARAMETERS 

In this section, the effect of the support stiffness and damping is 
treated through a boundary condition on the shear resultant. This 
resultant is assumed to act at the mean radius of the conical washer, 
implying that the contact area of the rubber is small in comparison 
to the area of the ceramic (see Fig. 3). Then, the forced harmonic 
response of the plate can be found from the solution to 



V 4 w — k 4 w = I 



-§«"", 0<r<r, 
0, r t < r < r c 



(12) 



where the wave number, k, is defined by 



* = P ~ (13) 

Due to the assumed cylindrical symmetry of the pressure, only 
axisymmetric solutions of (12) are sought. The plate is then divided 
into three regions: (i) < r < r, ; (it) r, < r < r, ; and (in) r, < r < r c . 
Solutions over these regions are pieced together by satisfying continuity 
(boundary) conditions on the transverse displacement, the slope, the 
radial bending moment, and the shear force at the radii r, and r, ; 
in addition, the homogeneous boundary conditions on radial bending 
moment and shear at the free outer edge of the plate are satisfied. 
At the effective support radius, 

Q(r.-) - Q(r.+) = WM. - fe.(«)] W (r.); (14) 

i.e., the net shear is opposed by a complex impedance that is propor- 
tional to the transverse displacement at that point. The impedance 
is composed of an inertia term, represented by co 2 M, , and a complex 
stiffness, written as a generalized Maxwell model 12 in the form 

fc.(w) = h'M + ik'/(u), (15) 

where 

k'.(u>) = k R + E r~^4-2 , (16a) 

wm = t r^h , deb) 

and k R , k n , and t„ are the equilibrium (rubbery) stiffness, an incre- 
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mental stiffness, and a relaxation time associated with an incremental 
stiffness, respectively. 

The technique that is used to solve this system, subject to the above 
stated boundary conditions, does not require the determination of the 
eigenfunctions of the problem. This procedure can be avoided since a 
particular solution is known: namely, 



k 4 D ' 



< r < r, 

(17) 
. 0, r t < r < r e . 

To this solution it is necessary only to add properly weighted homo- 
geneous solutions to satisfy the boundary conditions — in this case, 
ordinary and modified Bessel functions of zero order. The matrix 
inversion required to find the proper weights is carried out on a digital 
computer and the solution to (12) can then be determined as a function 
of the acoustic driving frequency. The voltage output is then found 
from the expression derived in the Appendix. 

Two numerical examples are solved in order to illustrate the pro- 
cedure and to compare the exact (lumped parameter) results with 
experiment. The transducer design is identical in both cases; the only 
difference is the conical washer material — in the first case, the PBD/SBR 
blend; in the second case, the PBD/CBT blend. Effective mass and 
stiffness are computed by procedures that were described previously. 
The comparison to experiment for the PBD/SBR blend is shown in 
Fig. 7 and a similar comparison for the PBD/CBT blend is shown 
in Fig. 8. 

The response comparison for both examples is favorable up to fre- 
quencies slightly above the first resonance; then, in both cases, an 
intermediate response peak is not captured by this model and the re- 
sponse peak at the next resonance is predicted to be much lower than 
shown by experiment. The location of this latter peak is, however, 
quite favorable. Perhaps the most disturbing feature of the comparison 
is the encouraging proximity of the analytical results to the design 
goal (see Fig. 1) — encouragement that is not borne out by the actual 
transducer performance. It seems apparent that other deformation 
mechanisms, not represented adequately by the lumped mechanical 
model of the conical washer support, are dominating the response at 
the higher frequencies. For this reason, a more exact model of the 
support structure is in order. 
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Fig. 7 — Voltage output vs frequency, rubber: 50 percent cis-4 PBD, 50 percent SBR. 



cr 

m -70 

CO 



EXPERIMENTAL 



— O— ANALYTICAL 







I I I I I I 



200 



300 400 



600 1000 2000 

FREQUENCY IN HERTZ 



4000 6000 10,000 



Fig. 8 — Voltage output vs frequency, rubber: 75 percent cis-4 PBD, 25 percent CBT. 
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V. FINITE ELEMENT ANALYSIS 

In order to understand the limitations of the lumped mechanical 
model of the rubber support, a finite element code 13 was exercised. 
The code was designed to dynamically analyze axisymmetric elastic 
solids subjected to arbitrary time-dependent loads and includes, as an 
option, the frequency and mode shape calculations for the solid. For 
this application, the ceramic was discretized into eighteen plate bending 
elements and the conical washers were discretized into three successively 
finer grids, with the most dense grid containing 92 quadrilateral con- 
tinuum elements. All materials were treated as elastic — the absolute 
value of the complex extensional modulus of the rubber was used — and 
Poisson's ratio was chosen to be either 0.45 or 0.49. 

A modal survey was then conducted for varying values of rubber 
extensional modulus. The four lowest resonant frequencies and their 
corresponding mode shapes were calculated for each modulus value. 
Typical results are shown in Figs. 9a-9d. These figures portray the 
influence of the rotatory inertia of the cantilevered frustum, which 
vibrates either in-phase or out-of-phase with the outer edge rotation 
of the ceramic. Note that the out-of-phase modes, Figs. 9a and 9c, 
are not strongly piezoelectrically active, whereas the in-phase modes, 
Figures 9b and 9d indicate substantial edge rotation with reference to 
central deflection. 

A composite plot of all the results obtained from the modal survey 
is shown in Fig. 10. This plot correlates well with the experimental 
results of Figs. 7 and S. Note that the results are only slightly de- 
pendent on the value of Poisson's ratio and on the discretization. 

VI. CONCLUSIONS 

With the knowledge gained from these three phases of analysis, a 
coherent set of design conclusions can be drawn. These recommenda- 
tions fall into two categories: (i) rubber material selection and (tt) 
conical washer design modification. In Figs. 11 and 12 the response 
variation of the ceramic transducer and its rubber supports is shown, 
as a function of environmental temperature, for the two different 
rubbers. 1415 Clearly, the rubber modulus is increasing too rapidly and 
the loss tangent is not holding an adequate value at the lower tem- 
peratures. In addition, the rotatory inertia of the unconstrained rubber 
is creating intolerable amplitude levels at the higher frequencies. In a 
recent investigation, 10 block copolymers cast from different solvents 
seemed to yield dynamic mechanical properties with desirable damping 
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Fig. 9a— Resonant mode shape, /o = 2250 Hertz, E R = 40 X 10« dynes/cm 2 . 
Fig. 9b— Resonant mode shape, /o = 2700 Hertz, E R = 40 X 10 s dynes/cm 2 . 
Fig. 9c— Resonant mode shape,/, = 5400 Hertz, E R = 60 X 10 6 dynes/cm 2 . 
Fig. 9d— Resonant mode shape, / 2 = S500 Hertz, E R = 80 X 10 6 dynes/cm 2 . 
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Fig. 10 — Finite element frequencies vs rubber extensional modulus. 
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Fig. 11 — Frequency response vs temperature change, rubber: 50 percent PBD, 50 
percent SBR. 
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Fig. 12 — Frequency response vs temperature change, rubber: 75 percent PBD, 25 
percent CBT. 



characteristics and relatively stable stiffness. Because of the dual 
transition (see Fig. 13), a styrene-butadiene-styrene block copolymer, 
obtained from solutions in carbon tetrachloride (c), toluene (T), ethyl 
acetate (E), and methyl ethyl ketone (M), has a sufficiently high loss 
tangent over a 200°C temperature range and also has a relatively 
constant modulus over a 130°C range. If the modulus is too high over 
this range, the washer design can be modified — thinner and taller 
cross section — to achieve nominal stiffness. A material tailoring program 
might produce a rubber which will help the transducer meet the design 
template. 

In addition to the improvement of rubber mechanical properties, 
the washer design should be altered in order to decrease, substantially, 
the rotatory inertia of the unrestrained cantilever section. One pos- 
sibility is shown in Fig. 14. An inverted vee-shaped design is depicted 
for the bottom washer that has three salient features: (i) a thinner 
cross section in order to maintain the modulus/stiffness ratio; (it) re- 
straint of the cantilevered section; and (Hi) a seating lip to aid in 
fabrication of the transducer. 
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Fig. 13 — E' and tan 6 vs temperature. 
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APPENDIX 

Voltage Output Calculation 

The voltage output of the bimorph ceramic microphone is calculated 
from the strain field under the following assumptions. 

The significant contributions to the component of the electric field, 
E z (r), perpendicular to the midplane of the plate are generated by the 
two strain components e rr and e 6 e ■ This assumes that the electric dis- 
placement field is negligible. Hence 

Z-E(r) = E z {r) = Z-(-h-S + ?• D) 

= /ll2£l2 "13 c 13 ) 

where Z is the unit normal vector to the plate, h the piezoelectric tensor 
relating strain, S, to electric field, and (5 the electric displacement 
field, D, to electric field. 17 However, for a ceramic poled in the Z direc- 
tion, h 13 = /f.,2 ; hence, 

E z (r) = -h l3 (e l2 + e 13 ) /jgv 

= —h l3 (t r r + €99). 

From thin-plate theory, the strain term is given in terms of the 
midplane displacement, iv(r), by the expression 

- + ---«&!('£)]: <*» 

thus 

The average electric field over the plated area is given by 

E nvP = \ f V ['' E z (r)rdrdO, (22) 

7T7"„ Jo J n 

where r„ is the electroded radius. 
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Substitution of (21) in (22) yields 

r v \ dr / r=rp 

The potential difference across one plate is therefore given by the 
line integral j h e E ave dz, but, since the two plates are series connected, 
the total voltage is given by 

n»-^*(fcL.. (24 > 

where h e is the thickness of one ceramic disk, and h = 2h e the total 
thickness. 
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